Large Scale Structure Formation of Normal Branch in DGP Brane World Model 
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In this paper, we study the large scale structure formation of the normal branch in DGP model 
(Dvail, Gabadadze and Porrati brane world model) by applying the scaling method developed by 
Sawicki, Song and Hu for solving the coupled perturbed equations of motion of on-brane and off- 
brane. There is detectable departure of perturbed gravitational potential from LCDM even at 
the minimal deviation of the effective equation of state w e a below —1. The modified perturbed 
gravitational potential weakens the integrated Sachs-Wolfe effect which is strengthened in the self- 
accelerating branch DGP model. Additionally, we discuss the validity of the scaling solution in the 
de Sitter limit at late times. 



PACS numbers: 



I. INTRODUCTION 



Dvali, Gabadadze and Porrati proposed a brane world 
model embedded in the Minkowski 5D dimension 
which has been taken as an alternative theory to explain 
cosmic acceleration. There are two distinct branches of 
this model: one is the self-accelerating branch (hereafter 
sDGP) and the other is the normal branch (hereafter 
nDGP). 

Deffayet discovered the accelerated phase in sDGP 
without cosmological constant [2 , which has been tested 
geometrically by many authors !, 0, H, IE 0- The de- 
tectability will be enhance d by the structure formation 
test 

The distinct evolution of perturbed gravitational poten- 
tial of sDGP has been studied by many authors at scales 
in quasi-static limit [H, [24], [25l. l2rj| and scales at nearly 
horizon size . sDGP strengthens the integrated Sach- 
Wolfe (hereafter ISW) effect, and leaves a detectable sig- 
nature on CMB power spectra and ISW-galaxy cross- 
correlation Q- 

There is no cosmic acceleration induced spontaneously 
in nDGP without introducing cosmological constant. 
But it's attracting features are that nDGP has a spec- 
trum leading to the stability in quantum level (27j . and 
that the effective equation of state w e g crosses w e s = — 1 
without phantom dark energy. nDGP has been tested 
by many authors geometrically as a natural theory cross- 
ing w e ff = — 1 line below [28j, |23, [30( . But no detailed 
study has been done to show how to formulate the large 
scale structure formation of nDGP in order to improve 
the detectability. 

We study the large scale structure formation of nDGP 
by using the scaling method developed by Sawicki, Song 
and Hu (26|. On-brane equations of motion are not in 
the closed form without the information of the gradient 
crossing the bulk direction. Off-brane equation of mo- 
tion should be coupled in order to supply this gradient 



of perturbations at the location of brane. We assume 
the scaling ansatz for the perturbation on the brane, 
and solve it iteratively until it converges. We find that 
the quasi-static limit is valid at scales relevant to ISW- 
galaxy cross-correlations. The modified perturbations 
weaken the ISW effect, which is the opposite behavior 
of perturbations found in sDGP. While we are preparing 
this paper, Cardoso, Koyama, Seahra and Silva indepen- 
dently derives the identical solution by using different 
approach [3l|. They verify that the scaling method is 
valid for both branches of DGP. 



II. LINEAR PERTURBATION EQUATIONS OF 
NORMAL BRANCH DGP 

Thin 4D brane is embedded in the 5D gravitational 
theory in the DGP brane model. The gravitational in- 
teraction propagates through the extra dimension, while 
all other physical interactions are confined on the brane. 
Tension is added to the matter-radiation energy on the 
brane in order to accelerate the cosmic expansion in 
nDGP at late times. The Einstein action is written as, 

s - -d? / ^ fl,5> - 2? / d '^ A "> 

+ J d A x^g(L m + L a ) , (f ) 

where L m denotes Lagrangian for matter-radiation en- 
ergy, L a denotes Lagrangian for brane tension, and bulk 
remains empty. The gravitational interaction scales are 
n 2 for 5D and /i 2 for 4D. The 5D gravitational interaction 
scale k 2 is a free parameter to be tuned to generate the 
observed cosmic acceleration. The ratio between both 
gravitational scales is defined by r c = k 2 /2/i 2 , which de- 
termines the critical scale of transition between 4D and 
5D gravity. 

The modified Friedman equation in DGP is given by 
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where p m denotes the matter density and p„ denotes 
the brane tension. There are two alternative choices 
for the sign convention of e, e = ±1. The choice of + 
sign is called as a self-acceleration branch which self- 
accelerates the cosmic expansion without cosmological 
constant. The other choice of — sign is called as a nor- 
mal branch. There is no self-acceleration in this branch 
without cosmological constant. In order to fuel cosmic 
acceleration, it should be aided by a term like the brane 
tension working as cosmological constant. 

Although it is less attractive to take nDGP as an alter- 
native candidate replaces LCDM, there is an interesting 
feature in this theory. The effective equation of the state 
of effective dark energy w e s is less than —1. The formu- 
lation is given by [3(3] , 



where A m and A^ denote the comoving perturbed en- 
ergy density of matter and Weyl fluid respectively (A m = 
S m — 3Hq m and Ae — Se — SHqe)- Signatures of terms 
with r c are reversed in nDGP. The geometrical anisotropy 
stress is generated by Weyl fluid anisotropy fluctuations, 
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Again, the overall sign is flipped in nDGP anisotropy 
stress equation. 

The conservation equations of matter fluctuations are 
given by 



T q m = -3$ 



(8) 
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where D, rc — l/4H 2 r 2 , and the expansion history factor 
E(a) = H(a)/H is 



E(a) = ^flm/a 3 + fl a + fl r 



(4) 



5D gravitational effect screens cosmological constant, 
which leads nDGP crossing w e e = — 1 line. 

The current constraint on the lower bound of w e ff is as 
low as w e g ~ —1.3 at a 95% confidence level by the com- 
bined test of Supernovae and WMAP [H, H3 ■ We vary 
w e g from —1.03 to —1.12. The cosmological parameters 
we use are: (u b = 0.025, w m = 0.128,^0 = 72, 10 10 A S = 
21.2, n s = 0.95, r = 0.09) from WMAP best fit [33]. The 
energy density of the brane tension fl a is varied with the 
demanding value of w g. 

We write the equations of motion of nDGP. There are 
two sets of equations. One is on-brane equations which 
are the projected Einstein equations and the conserved 
equations on the 4D brane. The other is off-brane equa- 
tion which is the propagation equation of perturbations 
through the bulk. 



q m = -* • (9) 
The conservation equations of Weyl fluid are given by 



k 2 

5e + H5e ^Qe = 

a 2 



where the source term is, 
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A complete set of on-brane equations of motion is given 
by combining all of these equations. Unfortunately, these 
on-brane equations of motion are not in close form, since 
the anisotropy equation is unknown due to the undeter- 
mined 7Te on-brane. We will show how to derive tte from 
the off-brane equation of motion which is called the mas- 
ter equation. 



A. On-brane perturbed equations of motion 



B. Off-brane perturbed equation of motion 



The perturbed metric on 4D brane is given by the New- 
tonian gauge, 



ds 2 = -(1 + 2^)dt 2 + a 2 (l + 2$)dx 2 



(5) 



The projected Einstein equations on 4D brane and the 
conserved equations give the on-brane equations of mo- 
tion. The projected Einstein equations provide a cou- 
ple of constraint equations; the Poisson equation and the 
anisotropy stress equation. There are two conservation 
equations of matter and Weyl fluid. Those equations 
complete the on-brane equations of motion. 
The Poisson equation is given by, 
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(6) 



In the maximal symmetry of five-dimensional metric, 
scalar perturbations can be expressed by a single scalar 
field called as master variable f2. The propagation equa- 
tion of fl through the bulk is given by [3J] 



O 



0. 



(13) 



where the 5D metric is d 2 s = —n(y, t) 2 dt 2 + b(y, t) 2 dx. 2 + 
dy 2 . The components of the metric of nDGP are given 
by 



n(y,t) = l-\^+H 
b(y,t) = a(l-H\y\) 



(14) 
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It is clear to see how the off-brane perturbed equation 
of motion helps us in solving the on-brane equations of 
motion by expressing Weyl fluid fluctuations in terms of 

n, 



1 
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where ' denotes the derivative in terms of y. The missing 
information on the brane is the gradient of Q crossing the 
bulk in tte equation. This gradient is supplied by solving 
Eq. 1131 the propagation equation of perturbations along 
the bulk. The coefficient of f2' term has the hidden e 
which should be reversed in nDGP. 

The perturbed Weyl fluid conserved equation can be 
re-expressed by using Eq. [15] in terms of f2, 

VL - 3HF(H)Q 



v 'a 2 K(H)r c r. 



2a 3 



k 2 K(H) 



2 A 

P Prn <-^rn 3 



where R denotes the gradient crossing the bulk at the 
location of brane, 



1 90 

R= midy' ly = ' 

and the coefficients F(H) and K(H) are 



(17) 
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(18) 



We complete our formulations for nDGP by switching all 
signatures in the terms including e. 



III. SCALING SOLUTION 

The complete set of equations of motion can be solved 
by using scaling method developed by Sawicki, Song and 
Hu [26[. The evolution of il is parametrized by the ex- 
ponent of a based upon the assumption of scaling ansatz 
on the brane, which enable us to intercommunicate be- 
tween on-brane and off-brane equations. We solve these 
equations iteratively until we get the converged pertur- 
bations. 

The master variable is parametrized by the exponent 
p on the brane as 



y=o= 



A(p)a p . 



(19) 



To take this as an initial condition, we can set the sep- 
arable function G evolving through the bulk coordinate 
y satisfying G \ y =o= 1. It is convenient to define a new 
variable x = yH, then the 5-dimensional scaling ansatz 
can be separable as, 



n(a,x) = A(p)a p G{x) . 



(20) 



The causal horizon of the propagation of perturbations 
through bulk is given by 



da' 



o a' 2 H(a') 2 



(21) 



The second boundary condition can be imposed at x — £ 
as G{x = £) = 0, until it hits the horizon in Gaussian 
normal coordinate y = 1/H (x = 1). We switch the 
boundary condition at x = 1 at £ > 1. It is fair to set 
this boundary condition since there are no other sources 
to propagate fluctuations beyond those horizons. 

We put this scaling ansatz into the off-brane equation 

Eq.rjn 



d 2 C AC 

+ B{x)— + C(x)G = , 
ax z ax 



(22) 



(16) where the coefficients are given by 



A(x) = {l-x)[l-x{\ + 2h)) 

B{x) = -2xhp - 2xh 2 - xti - (1 + h) + x(l + h) 2 



x 2 h 



h' + h 2 + h 
l-x(h + l) 
, [1 - x(l + h)] 2 



3xh 



[1 - x(l + h)} 
1-x 



1-x 



C(x) 



h' + h 2 + h 



3p 



-p -hp-xp- 

1 — x(h + 1) 

l-x-xh [l~x(l + h)} 2 k 2 



l~x 



{l-xf 



a 2 H 2 



(23) 



where h denotes (d-ff/d In a)/_ff . When we solve Eq. 
it gives the gradient coefficient R in terms of G, 



R 



dG I 

dx \ x=0 

G L=o 



(24) 



Then we can solve the on-brane equation Eq. 1161 

We start the iteration from solving the off-brane equa- 
tion with the trial value p. We use p = 4 which is the 
expected value of p during matter domination epoch. It 
returns the value of R, and now the on-brane equation 
can be solved. The on-brane equation gives more realistic 
value of p, 



P 



dlnft 
dlna 



(25) 



Then the same process is repeated until the solution con- 
verges. 
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FIG. 1: The perturbed potentials ($ — ^E')/2 in quasi-static limit. 
The effective equation of state w e ff varies with f2 CT with fixing all 
other cosmological parameters. The dotted curve represents (<E> — 
9)/2 of LCDM, and the solid curves represent (<t> - \P)/2 with 
varying w e f[ from —1.03 to —1.12. 



FIG. 2: The perturbed potentials (<t> — \&)/2 at scales of nearly 
horizon size are plotted. The solid curves represent the quasi static 
solutions, and the dash curves represent the dynamic scaling solu- 
tions scales at k = 10 _3 Mpc _1 and k = 5 X 10 _4 Mpc~ 1 . 



In the quasi-static limit, the spatial gradient term dom- 
inates over all other quasi-statically varying terms. Then 
tte is simply related to 8e by constraint equation to close 
on-brane equations of motion, 

3 

The perturbed potentials can be written in the simple 
form, 

k ^=4. G (l-±) P A m (27) 



^ = -4.G(l + -) P A m , (28) 

where 

(3=1 + 2Hr c (l + . (29) 

The reversed sign in (3 leads to (1 + 1/3/3) > 1. Then the 
Newtonian potential well becomes deeper in nDGP than 
in LCDM, and matter fluctuations falling into that well 
are enhanced. It is the opposite behavior of perturbations 
observed in sDGP where matter fluctuations decay due 
to the shallower Newtonian potential well. The dynamic 
solution from iteration routine precisely reproduces this 
result. The quasi-static limit is valid for nDGP. We plot 
the solution in Fig. [T] with varying w e s from —1.03 to 
-1.12. 



For the modes in crossing the horizon, the quasi-static 
approximation is invalid. Weyl fluid anisotropy is in- 
fluenced by all other terms including the bulk gradient. 
The dynamic solution starts to depart from the solu- 
tion of the quasi-static approximation. Its behavior is 
presented in Fig. [2] at larger scales of k = 10 _3 Mpc _1 
and k = 5 x 10 -4 Mpc . The dynamic solution departs 
from the quasi-static solution at scales of k smaller than 
k < 10 _2 Mpc _1 . The ISW effect on CMB power spectra 
is influenced by this departure, while the ISW effect on 
ISW-galaxy cross-correlations is well described by using 
quasi-static approximation only. 

We test the numerical stability by using Bcrtschinger's 
equation. Bertschinger derived the dynamic equation for 
the perturbed potential at super-horizon modes based 
upon just two conditions, the Friedman universe and the 
conservation of energy. As long as these two conditions 
are not broken in the theory, the following equation is 
always valid at the limit of k — > [H, [36| , 



^_^_|_$'+^_| r j$ = (30) 



here ' denotes the derivative in terms of In a. We take 
the numerical anisotropy stress in the limit of k — > and 
feed it into Eq. 1301 The numerical perturbed potentials 
from the scaling method are consistent with the solution 
from Eq. H3 
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IV. THE VALIDITY OF SCALING METHOD IN 
ASYMPTOTIC LIMIT 



The growing solution of Q with this p leads to the growing 
anisotropy. With the choice of c ~ 2, $ + becomes 



We derive the stability condition for nDGP in the 
asymptotic limit. In the future, nDGP approaches to 
the de Sitter universe. All modes are expelled outside 
the horizon, £ becomes infinity, and k/(aH) 2 <C 1 limit is 
available for all modes. The critical distance r c is greater 
than the constant asymptotic particle horizon l/iJ(oo), 
which makes c = H(oo)r c > 1. In the de Sitter limit, the 
master equation is written as 



n E . i 

+ oc — ~ A + a 
a 



(37) 



d 2 n 

dy 2 



2H dn 



H 2 



1 — Hy dy (1 — Hy) 2 \d\na 2 din a 



d 2 n 



dn 



0(31) 



Two boundary conditions, Q(y = 0) = 1 and Q(y = 
1/H) = 0, leads to the solutions trivially in the following 
forms, 

Q = A(p)a p [B+(l- Hy)*+ +B-(1 — Hy)*- ] . (32) 

Since (1 — Hy) is zero at y = 1/H, the negative exponent 
of (1 — Hy) will violate the normalizability. Thus only 
the positive exponent will be acceptable with B + = 1 to 
satisfy the other boundary condition at y = 0. 

The solution of master equation, Eq. [331 gives us the 
bulk gradient at y = 0, R = —q+. Then the conservation 
equation of master variable on the brane in de Sitter limit 
is simply written as, 



9 1 1 

q+ + -q+ - - = o , 

c c 



(33) 



where we ignores the contribution of matter perturba- 
tions which converges to constant at late times. The 
positive solution q + is 



1 1/1 4 



(34) 



It gives the evolution equation of on the brane at late 
times, 



P - 3p 



"1 








c 





= 0. 



The growing mode has the exponent p as, 




4g 4 



(35) 



(36) 



But it is a misleading conclusion of asymptotic behav- 
ior of nDGP by using scaling method improperly at late 
times. The negative R violates the renormalizability con- 
dition [37:]. It indicates that the real exponent of q can 
not be a solution, and it should be imaginary. Then 
we are not able to ignore B_ by regularity condition at 
y = 1/H. Since the boundary condition at y = 1/H 
does not determine any coefficient of B's, the solution 
of nDGP at asymptotic region is unbounded. Thus the 
scaling solution is invalidated in the de Sitter limit [38| . 
We can use the scaling method when the matter pertur- 
bations keep growing and attract fl to grow, up to the 
present time or near future until hitting de Sitter limit. 



V. CONCLUSION 

We have solved the dynamic equations of nDGP by 
using the scaling method. A detectable departure from 
LCDM is found with w e g below — 1. The matter fluc- 
tuations are enhanced due to the deeper potential well 
in nDGP than in LCDM, which weaken the ISW ef- 
fect. We are able to probe nDGP by detecting this 
weakened ISW effect by ISW-galaxy or ISW-weak lens- 
ing cross-correlation. It is highly recommended to select 
high redshift bins where the modified perturbed poten- 
tial increases to provide the negative cross-correlation. If 
found, this would be strong evidence for the presence of 
the nDGP model [H. 
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